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Abstract. - The helix-coil transition is modified by grafting to a surface. This modification 
is studied for short peptides capable of forming Q-helices. Three factors are involved: (i) the 
grafting can induced change of the boundary free energy of the helical domain (ii) the van 
der Waals attraction between the helices and (iii) the crowding induced stretching of the coils. 
As a result the helix-coil transition acquires "all or nothing" characteristics. In addition the 
transition temperature is elevated and the transition itself sharpens as the grafting density 
increases. 



During the past two decades the physics of polymer brushes formed by terminally anchored 
chains were studied extensively |l], |[ |j. Most of the research effort dealt with brushes of 
flexible, synthetic polymers, devoid of internal degrees of freedom. In contrast, this letter 
concerns brushes formed by biopolymers capable of undergoing a cooperative helix-coil tran- 
sition. It is motivated by two experimental observations. First, the promotion of the adhesion 
and spreading of cells by brushes of collagen model peptides Q. Second, membrane fusion 
induced by model fusogenic peptides grafted to vesicles Q. In both cases, the function of the 
short peptide chains was correlated with a helical state. Furthermore, the helical state was 
favored by the grafting. With this in mind we present a highly simplified theoretical model 
for the helix-coil transition in brushes of short, laterally immobile peptides. In particular, 
we discuss the transition temperature T t and the width of the transition AT as a function 
of grafting density, S^ 1 . We focus on the simplest situation, of short neutral homopeptides 
forming a single stranded a-helix. As we shall see, the grafting of the chains can lead to 
qualitative modifications of the helix-coil transition: (i) in marked distinction to the case of a 
free peptide, the helix-coil transition in an isolated grafted peptide can acquire "all or nothing" 
characteristics; (ii) T t is elevated with the grafting density while AT decreases; (iii) eventually, 
for high grafting density the transition may take place as a first-order phase transition. These 
distinctive features arise because of three factors: The lower configurational entropy of the 
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monomer at the grafting site favors helix formation even in isolated chains. In a brush, the 
helical state is also promoted by the crowding induced stretching and by the attractive van 
der Waals interaction between the helices. 

In a a-helix |6[ 0, || the ith monomer forms i?-bonds with the (i — 3)th and the (i + 3)th 
monomers. Overall, a helical domain consisting of n monomers contains n — 2 ii-bonds. It is 
convenient to consider the chain as a sequence of bonds. Choosing the coil state as a reference 
state, each helical bond is associated with an excess free energy Af reflecting the formation 
of ii-bonds, the accompanying change in solvation and the loss of configurational entropy. 
Af is a function of the temperature, T. A/ vanishes at the transition temperature T*, while 
Af < when T < T* and A/ > for T > T*. The terminal bonds, at the boundary of 
the helical domains, are associated with an additional free energy penalty A/ t . This arises 
because the terminal bonds lose their configurational entropy but do not contribute ii-bonds. 
Aft plays the role of an interfacial free energy associated with the helix-coil boundary. It is 
customary to formulate the theory of the helix-coil transition in terms of the Bragg- Zimm 
parameters s — exp(— Af/kT) and a — exp(— 2Af t /kT). a is independent of T and is 
typically of order of 10~ 3 — 10~ 4 , depending on the identity of the amino acid residues forming 
the peptide. On the other hand, s is a function of T and in the vicinity of T* it varies 
as s — 1 ~ (T* — T)/T*. In an infinite chain, the plot of the fraction of helical bonds, 9 
vs. s is sigmoid and the width of the transition is AT ~ T^er 1 / 2 . Since the chain is one 
dimensional object, a first order phase transition is impossible. As a result the chain consists 
of an alternating sequence of helical and coil regions. The minimal size of a domain comprising 
a helical and a coil regions, as obtained at T», is er -1 ' 2 . When the polymerization degree, N is 
much larger than tr -1 / 2 the chain may be considered as infinite. It consists of a large number 
of domains and the width of the transition, a 1 / 2 , is independent of N. On the other hand, 
when N < cr -1 / 2 the chain incorporates typically only one helical region and the width of the 
transition is of order 1/N. The statistical physics of such short oligopeptides are well described 
by the "one sequence approximation" where the chain is assumed to contain a single helical 
region Q . The customary formulation of this approximation utilizes the appropriate partition 
function. For the purposes of our discussion it is convenient to utilize the corresponding free 
energy. The chemical potential of a chain supporting a helical domain consisting of n bonds is 
/j,q (n) — nAf + 2 Aft — kT In (AT — n — 1) where the last term is the entropy associated with the 
placement of the helical segment along the chain. Altogether there are N — 2 bonds and the 
first monomer of the helical segment of length n can be placed in any of the (N — 2) — n + 1 
sites. For n = we have Ho{0) — 0. The fraction of chains with a given n, p n , is determined 
by minimizing £1 = Do + M^2 n p n [^{n) + kT\np n } subject to the constraint ^2 n p n = 1. This 
leads to /i = fj,o(n)+kT 'lnp„ = const' since all species coexist in equilibrium. For our discussion 
it is sufficient to consider the dominant term in f2, as specified by d[io(n)/dn = 0. When 
Af > the minimal fio(n) is /io(0) = and the majority of chains contain no helical region. 
On the other hand, for Af < the probability distribution peaks at = — 1— | kT j Af j. 
Accordingly, n* attains its maximum value, n* = N — 2, when | Af \= kT. The discussion 
as presented above applies to free chains. Two features are of special importance for future 
reference: (i) Af t is independent of the position of the helical region and, as a result (ii) all 
placements of the helical segments are equally probable. Because of the associated entropy n» 
is lower than N - 2 for > Af > -kT. 

The situation described above is modified significantly when the chain is grafted, terminally 
anchored, to a surface. We first discuss the case of a single chain grafted to a flat surface. 
The grafting gives rise to two effects. First, the overall number of configurations of the 
chain is reduced because of the presence of an impenetrable wall. As a result, the overall 
configurational entropy is reduced by AS ~ In JVfl . This is expected to lead to a small changes 
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in s and a which we will ignore since they do not give rise to qualitative effects. Second, the 
two boundaries of the helical region are no longer equivalent because of the presence of the 
surface. The surface effect is expected to decay with the length of the coil region separating 
the helical region from the wall. The maximal effect is attained when the helical region is 
initiated at the wall. The penalty of the "free" boundary remains Aft. However, the terminal 
penalty at the wall, Af w , can be different. The difference can arise from two sources. First a 
modification of the torsional potential for the terminal monomer and the consequent reduction 
of its configurational entropy. Second, change in the electrostatic interactions experienced by 
the terminal monomer at the grafting site. As a result, the Bragg-Zimm parameters are 
modified and for a helical sequence at the wall a g — exp[— (Af w + Af t )/kT]. In the following 
we consider the physically plausible case of Af w < Af t . For simplicity we assume that Af w 
is attained only when the helical sequence begins at the grafting site. In this situation a 
helical sequence of length n may assume two states: (i) When the terminal monomer is not 
at the wall, the helical sequence can be freely placed at the available remaining sites. The 
"free" sequence is associated with /i/o = nAf + 2A/t — kT\n(N — n — 2). (ii) A "bound" 
state when the terminal monomer is at the grafting site and fibo = nAf + Aft + Af w . The 
two states are equally stable when fifo = fibo- This condition is satisfied for n-f given by 
n-f = N — 2 — exp[(A/t — Af w )/kT]. For n > n-f the bound state is favored. Accordingly, the 
bound state is always dominant when n+ = 1 or ln(A — 3) < (Aft — Af w )/kT. Furthermore, 
among the bound states, the fully helical state is always the most stable i.e., Hbo(N) < Hbo( n ) 
whenever Af < 0. For this choice of N the coil state remains the most probable state, n» = 
when Af > while for any Af < the most probable state is fully helical, n* = N — 2. Thus, 
the grafting of short chains modifies the nature of the helix-coil transition that acquires "all 
or nothing" characteristics. In the remainder of this letter we focus on brushes of such chains. 

Our earlier considerations concerned a single grafted chain. With them in mind we model 
a grafted layer consisting of many chains as a mixture of chains in a coil state and chains in 
a fully helical state. In this case, additional contributions come into play. One is the van der 
Waals attraction between the rods. Another important contribution is the crowding induced 
extension of the coils in the brush regime. To quantify this picture it is necessary to introduce 
further assumptions concerning the orientation imposed on the helices by the grafting. For 
simplicity we will consider the case of grafting sites enforcing perpendicular orientation of 
the helices with respect to the surface. We further limit the discussion to the case of rigid 
junctions that is, the helices cannot bend at the grafting sites. Thus far, our considerations 
were directed at chains with N < er~ 1//2 w 10 2 . At this point it is necessary to limit our 
discussion further to chains with TV < « (a/a/,.) 2 where a is the effective monomer size in 
the coil state while ah = 1.5 A. is the projection of a residue on the axis of the helix. While 
a depends on the identity of the residue, a/ah > 1 and typically A» > 50. This choice is 
necessary because the length of a fully helical chain, L, is smaller than the radius of the 
corresponding random coil, Rq rs A fl / 2 a when N < JV* while L > Rq if N > TV*. Thus, in 
one case the helices are submerged in the brush formed by the coils while in the other the 
helices can be partially exposed (Figure 1). This choice of the N range is dictated by the 
experimental systems motivating this work |^|, |^] . With the system fully specified we are in a 
position to write down the free energy per chain, F. For simplicity we use a modification of the 
Alexander model [[to] , assuming that the coils are uniformly stretched and that concentration 
profile of the monomers of the coils is step-like. We consider a planar grafted layer such 
that the area per chain is £ and the fraction of chains in a helical state is 9. F consists 
of four terms, F = Fq + 9F v dw + (1 — 9)F CC + 8Fh c - Fq reflects the contribution of the non 
interacting chains F « d(NAf + Af t + Af w ) + kT[9 In 8 + (1 - 9) ln(l - 9)] where the last term 
is the mixing entropy associated with the possible placements of chains in the two possible 
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Fig. 1. - A brush of short chains undergoing a helix-coil transition. The helical chains are short 
compared to the span of the coils and are thus fully embedded in the brush formed by the coils. 



states. The van der Waals attraction between two oriented helices is modeled as the interaction 
between two parallel rods of diameter d. When the distance between two neighboring rods, 
D, is small in comparison to L, their length, F v dw ~ — L/Z3 5 [pT| while in the opposite limit 
Fvdw ~ —L 2 /D 6 . In the first case, an element of the rod experiences, in effect, an interaction 
with an infinite neighboring rod. In the second limit, F vc iw reflects the sum of all pairwise 
interactions between the elements of the two rods. The van der Waals energy of a rod within 
the grafted layer, where the minimal distance between two neighboring helices is (AE/Trff) 1 / 2 , 
is 

ALd 4 fn9\ 5/2 Itt0L 2 

F v dw = — arctanv — — — (1) 

W 18 \AZJ V 4E y ' 

where A is the Hamaker constant. In a "normal" brush, the overlap between the coils gives 
rise to chain extension along the normal to the surface. The stretching reduces the number 
of monomer-monomer contacts at the price of an extension penalty. In our situation the 
picture is somewhat different. First, it is necessary to allow for both coil-coil and coil-helix 
interactions. Second, only the coils are extendible. Finally, for our choice of N, the helices are 
fully immersed in a brush of coils. F cc allows for the elastic free energy of the coils as well as for 
monomer-monomer interactions involving the coils. This contribution has the form of the free 
energy of a brush formed by coils that is, F cc /kT « H 2 /Na 2 + v cc N 2 a 3 (l — 9) /HI where H is 
the thickness of the brush and v cc is the virial coefficient for monomer-monomer interactions 
involving coils p2[. The first term reflects the Gaussian stretching penalty while the second 
allows for the repulsive monomer- monomer contacts. F cc must be supplemented by Fh c that 
reflects the interactions between the coils and the helices, Fh c /kT rs w/j C iV 2 a 3 (l — 9) /if£. 
Here Vhc is the virial coefficient for the binary interactions between helical monomers and coil 
monomers. For simplicity we consider the case of Vhc = v cc = v. In this case the equilibrium 
thickness of the brush H eq /a ~ N(ya? /H) 1 ' 3 , as specified by dF/dH — 0, is independent of 6 
so long as 9 is small enough to ensure overlap between the coils p^j]. For such a choice of 9 the 
equilibrium form of the free energy per chain, as obtained by substituting H eq into F, is 

"\f ejAft + Af w ) yy7 , fo 2 ^ 5 ~ 



NkT kT NkT \ S 



+l[01n0 + (l-0)ln(l-0)] (2) 
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where A ps wWa^' F v dw in @ is approximated by the ir9L 2 /£ 3> 1 limit of (|TJ). As a result, 
F v dw is overestimated for 9 <C 1. This does not affect our analysis since the contribution of the 
van der Waals attraction in this regime is negligible. The corresponding spinodal condition, 
d 2 F/d6 2 = 0, leads to f(9) = 9 5 / 2 (l - 6) ps (4/35)(E/a 2 ) 5 / 2 (AA0 _1 revealing a critical point 
at 9 C = 5/7 where f(6) exhibits a maximum. The critical grafting density is thus specified by 



E c /a 2 ps A 2/5 N 2 / 5 



(3) 



and for E < E c the helix-coil transition within the layer proceeds as a first-order phase 
transition. For lower grafting densities, £ > £ c , the helix-coil transition is cooperative but no 
phase transition is involved. In this last regime it is of interest to characterize the dependence 
of the transition temperature T and the width of the transition AT on E and N . To this end 
it is helpful to recast the equilibrium condition, dF/89 = 0, in the form 



1 - 



K(9) 



N 

s <T g exp 



-XN 



9a^ 
~E~ 



5/2 



N 



va 



2/3' 



(4) 



where K(9) is the equilibrium constant governing the ratio of helices and coils. A rough 
idea concerning the transition is obtained by identifying it with the condition 9 = 1/2 i.e., 
K(9) = 1. It is helpful to consider first the case of non interacting mushrooms, £ — > oo, when 

K{9) — s N <7 g . Since s in the vicinity of T* is given by s ps exp —r — 1\ where r is a 

phcnomcnological constant, this condition leads to 



T t (oo) ps T* 



In a „ 

Nr 9 



< Tl 



(5) 



Similarly we define the width of the transition as AT = T_ — T + where T + is the temperature 
for which 9 + = 9/10 while T_ is the temperature corresponding to 6>_ = 1/10 Jul. The ratio 



K{9+)/K(9-) = (s + /s_) JV = 81 may be rewritten as In 81 ps NrAT(oo)/T* or 
By using the same argument for a brush with E <C i?Q we obtain 



T = T(oo) + 



7, 



-1 

2E J 
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and 



AT = AT(oo) - - — A ( 

v ' 2 r V 10S / 



5/2 



(G) 



(7) 



(8) 



As E decreases, T t increases while AT decreases. In other words, the stability of the helical 
state grows and the transition becomes sharper when the grafting density increases. Eventually, 
at the critical point, AT = thus signaling the onset of a phase transition. Note that within 
our model the decrease of AT is due only to F v dw while the increase of T t results from 
contributions from F V( iw as well as the brush penalty. This is because our choice of v^ c — v cc 
leads to a 9 independent H eq . In the general case AT should reflect both contributions. 

The analysis presented above focused on the case of short, immobile, grafted chains that 
form a single-stranded helix. In this system T t increases while AT decreases as the grafting 
density grows. The effect of grafting on a, the van der Waals attraction and the crowding 
induced stretching should however play a role irrespective of the precise specifications of the 
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grafted layer. Additional factors may contribute when chain mobility, multiple stranded helices 
and longer chains are considered. Lateral chain mobility may give rise to in-plane phase 
separation driven by the van der Waals attraction. The discussion of multiple stranded helices 
should allow for the effect of loops on er||. It may also be necessary to allow for loss of 
translational entropy due to the formation of a multiple helix. The "all or nothing" model 
described above is clearly limited to brushes of short chains. A discussion of brushes formed 
by long helix-forming chains should allow for helix-coil coexistence on each of the polymers. 
While the discussion focused on brushes of peptides, it should be noted that similar situation 
is encountered in DNA chips undergoing hybridization jn|. Finally, it is of interest to note that 
our considerations are somewhat similar of the analysis of the coupling between the helix-coil 
transition and the onset of liquid crystalline order in peptides solutions |l6) and in a collapsing 
chain[[l7|. From this perspective, the distinguishing features of the brush case are due to the 
grafting modification of a and to the crowding induced chain stretching. 

The authors benefited from instructive discussions with M. Tirrell and D. Leckband. 



REFERENCES 
[1] S.T. Milner, Science, 251 (1991) 905. 

[2] A. Halperin, M. Tirrell and T.P. Lodge, Adv. Polym. Set. 100 (1992) 31. 

[3] I. Szleifer and M.A. Carignano, in Advances in Chemical Physics (Vol XCIV, Ed. I. Prigorgine 
and S.A. Rice) 1996. 

[4] (a) P. Brendt, G.B. Fields and M. Tirrell J. Am. Chem. Soc. 117 (1995) 9515, (b) Y.-C. Yu, P. 

Brendt, M. Tirrell and G.B. Fields J. Am. Chem. Soc. 118 (1996) 12520, (c) G.B. Fields et al 

Biopolymers (Peptides Set.) 47 (1998) 143. 
[5] E.-I. Pecheur et al Biochemistry 37 (1998) 2361. 

[6] T. M. Birshtein and O. B. Ptitsyn, Conformations of Macromolecules (John Wiley, New York) 
1966. 

[7] D. Poland and H. Scheraga, Theory of Helix-Coil Transitions in Biopolymers (Academic Press, 
New York) 1970. 

[8] A. Yu. Grosberg and A. R. Khokhlov, Statistical Physics of Macromolecules (AIP Press, New 
York) 1994. 

[9] E. Eisenriegler Polymers near Surfaces (World Scientific, Singapore) 1993. 
[10] S. Alexander J. Phys (France) 38 (1977) 977. 

[11] J.N. Israelachvili, Intramolecular and Surfaces Forces (2nd edn. Academic Press, London) 1991. 
[12] The van der Waals interactions between the monomers of the coils are allowed for in the Flory \ 

parameter or, equivalently, in the second virial coefficient v cc = (1 — 2\) as discussed in P.G. de 

Gennes Scaling Concepts in Polymer Physics (Cornell University Press, Ithaca) 1979. 
[13] In the vicinity of 9 — 1 the coils no longer overlap and some of the helices do not interact with 

the coils. Both F cc and Fhc assume different forms in this regime. 
[14] The definitions of T t and AT are clearly somewhat arbitrary. The rigorous definitions, in terms 

of the curvature of the 6 vs. T plots, do not allow for simple analytical solutions. 
[15] D. J. Graves Trends Biotechnol. 17 (1999) 127. 
[16] Y. H. Kim and P. Pincus Biopolymers 17 (1979) 2315. 
[17] A. Yu Grosberg Biophysics 29 (1984) 621. 



